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Abstract 

The equilibrium equations for wide, developable, Mobius strips undergoing large deformations 
have recently been derived, and solved, numerically. We use these results to compute the eigen- 
values and eigenstates of non-interacting electrons confined to Mobius strips of linking number up 
to 1.5 and of arbitrary width. The inverse participation ratio is used to show that electrons are 
increasingly localised to the higher curvature regions of the higher-width structures, where sharp 
creases could form channels for particle transport. Our geometric formulation could be used to 
study transport properties of Mobius strip and other components in nanoscale devices. 



I. INTRODUCTION 

The eigenstates of a particle confined to an inextensible Mobius strip have previously been 
calculated, assuming the Schwarz parametrisation [l| as an approximation to the centreline 
of the strip, for the single, large, aspect ratio, L/2w = 30, corresponding to the dimensions of 
ribbon-shaped crystals of NbSes Mobius shell structures that have recently been fabricated 
B. For this aspect ratio, usmg the Schwar. parametrisatiorr. the berrding errergy is easily 
minimised, after performing the Wunderlich reduction [3|, |J| to a one-dimensional integral 
over the centreline of the strip. The constraint that the length of the strip must equal 
L allows one of the three Schwarz parameters to be eliminated, giving a two-dimensional 
visualisation of the bending energy in the remaining two Schwarz parameters. For L/2w = 
30, one minimum is clearly visible. For smaller aspect ratios (larger widths), however, as 
considered in [3|], e.g. L/2w = 1, there are three candidate minima, none of which satisfy 
the constraint which ensures the generators do not intersect within the strip l|, |J], I^X"")! — 
w~^,u G [0, 27r] (for rj, see below), indicating the limitations of the Schwarz parametrisation 
even for modest strip widths. 

We note that recently |5| , wave functions have been calculated for graphene nanoribbons 
of varying Lk, again for a single, large, value of L/2w = 6.8. After a classical geometry 
optimisation for a discrete (atomistic) model, an annealing simulation searches for lowest 
energy structures of a classical dynamics. The electronic structure is simulated using a 
semi-empirical Hamiltonian. 

In order to enable, and generalise, the calculation of energy eigenvalues and surface eigen- 
states of a non-relativistic quantum particle to arbitrary widths, therefore, we use the formu- 
lation in 3[] which rigorously minimises the Wunderlich elastic bending energy functional ^ 
to obtain the exact shape of Mobius strips of linking numbers Lk = 0.5, 1.5. Consideration 
of the large width behaviour of the inverse participation ratio provides evidence for curva- 
ture trapping in these periodic structures normally seen in infinite domain Gj or disordered 
systems [7|. 



II. THE QUANTUM MECHANICS OF A PARTICLE BOUND TO A SURFACE 

In the formulation 8j for an infinitely strong squeezing potential which attaches a particle 
to a surface, the transformation of the wave function, ^z^, to x = v^V'; with F = 1 — 
2Mv? + K{u^Y, leads to a separation of the three-dimensional Schrodinger equation for 
X = Xt{u^,u'^)Xn{u^), where u^,u^ are the coordinates embedded in the surface, and u^ is 
the distance perpendicular to it, giving for the surface wave function Xt'- 

-^(A + M')xt = Ext, (1) 

where A is the Laplace-Beltrami operator on the surface. For a developable (inextensible) 
surface, the Gaussian curvature K is zero, and the mean curvature is easily calculated, using 
the coefficients of the first and second fundamental forms of the surface, as 

M = -^^±^. (2) 

21 + t7]' ^ ^ 

If r(s) is a parametrisation of the centreline of the strip, then 

^is,t) = ris)+t[his)+vis)t{s)] 

(3) 
t{s) = ri{s)K{s), s e [0,L],t e [—w,w] 

is a parametrisation of a strip of length L and width 2w, where t is the unit tangent vector, 
b the unit binormal and where k and r are respectively the curvature and torsion of the 
centreline [9i| which uniquely specify (up to Euclidean motions) the centreline of the strip as 
well as the Serret-Frenet basis vectors t, n, b. The surface is developable and is completely 
determined by the centreline of the structure. Developing the surface into a rectangle with 
rectangular coordinates (m^,m^), given by 

u^ = s + t7]{s), u'^ = t, (4) 

A in (pP) is then the usual Laplacian and M^ provides a quantum potential well. The 
boundary conditions are 

x(m\ m^ = -w) = x(m\ u"^ = w) = 0, 

(5) 
x{u' + 2L,u')=x{u\u'), 



where the latter is the requirement of the single- valuedness of the wave function [lOJ] 



known as the periodic (or Born-von Karnian) boundary condition 12(] 



xlO 



x10 



10 










1 






'A 


4 














0.5 






j\ 


2 








y 5 






f\ ^ 






/ '■ 


V 




^-^^ 


^ 






J 




\ 


-0.5 






-2 



















H 






A 










- Exact 


1 


2 1 


—4 

2 C 


) 1 


2 

x10"^ 




Sch 


warz 


s 


X 


10"° 




b 


X 


10"° 




s 



FIG. 1: Lk = 0.5,L/2u; = 30 



III. NUMERICAL RESULTS AND DISCUSSION 



A. Single twist Mobius strip: Lk = 0.5 

Shown in Fig. [1] are k, 77, rl as a function of the arc-length, s, in [3|] for the exact shape 
of a free-standing Mobius strip, characterised solely by its bending energy, for L/2w = 
30, compared with the Schwarz parametrisation. Clearly n and t] have an even reflection 
symmetry about s = L/2 whereas rj' has odd reflection symmetry which through (j4]) induce 
the following transformation in (m^,m^): 



u 



L-u\ 



2 2 

u —>■ —u . 



(6) 



The Hamiltonian H = — ^(V^ -|- M^) is invariant under this transformation since M^ is. 
Therefore non-degenerate states are such that 



X{u ,u ) = ±x{L-u ,-u ). 



Similarly Ti is invariant under the transformation 



u^ -^u^ + L, 



u 



-u 



itself induced by the transformation s ^ s + L, with corresponding parity eigenstates 



(7) 



(8) 



(9) 



Eq. ([9]) can be recognised as the Bloch or Floquet theorem for periodic potentials: 

x{u' + Ly)=e^^'^x{u\-u'), (10) 



which, given the periodic boundary condition, gives kL = nn, with n an integer. Eq. (fTOj) 
then gives ([9]). Thus four different symmetry eigenstates are considered: even and odd 
reflection symmetry in the hne u^ = L/2, and even and odd symmetry under translation 
by L, each with m^ —>■ —u^. Reflection symmetry, ([7]), allows the domain for the numerical 
computation to be reduced to half the strip. We used a finite-difference (FD) scheme, 
calculating the eigenvalues and eigenstates with MATLAB. Thus in the FD scheme, using 
^ for u^ < 0, M2(u\u2) = M2(2L - |m^|,m2) = M^{L - \u^,-u'^), (as b has changed 
sign and we are exactly under the surface from where we started), = M^(L + |m^|, —u^)-, (by 
symmetry), = M'^{\u^\^u'^) (as b has changed sign). This enables us to distinguish between 
particles on the strip from those in the strip (cf. [13.] for a flat Mobius lattice model), the 
former allowing negative parity eigenstates under translation by L. Without imposing (jTj), 
we find the FD scheme skips some eigenstates. The results were also checked with finite 
elements. The requirement v^ -^ —u"^ for invariance of M ^ under the s -^ s + L translation 
is because we use a continuously varying t, n, b frame moving along the centreline, changing 
to an anti-Frenet frame to avoid a Frenet frame flip where k = (n ^ — n), which therefore 
defines the co-ordinate system used in u^ and m^. In addition, we therefore require that 
7] -^ —7] and 7]' -^ —Vj and k -^ — k, under s ^ s + L, to define the t, n, b frame used. 
Otherwise, if the Frenet frame is used throughout, it flips under translation, and u^ -^ —u^ 
would not be required. 

Shown in Fig. [2] are contour plots of M^ for four increasing widths, taken from |3|, which 
form the potential wells which scatter the standing waves of electrons confined to the strip, 
i^or clarity, the outer boundary is omitted, as there are singularities there for |?7'(s)| = Xjw 
3|]. As w increases, creases are formed in the Mobius structure, which a quantum particle 
experiences as deepening potential wells which lower its energy. For low enough aspect 
ratio [Ljlw = lOvr/lS), negative energy eigenstates appear, the usual signature for bound 
states. Shown in Fig. [3] (left) is the dimensionless energy E' = {h^w'^ /2m)E versus w. By 
comparison with the non-geometric (flat) strip (obtained with M^ in ([T]) set to zero), we 
expect the energy to decrease with increasing width. There are therefore two competing 
effects. Multiplying E by w'^ makes E' increase with width initially, but for larger widths, a 
decreasing E' means that E decreases faster than w"^. The figure follows the lowest energy 
eigenstate for each parity with increasing width, with either one or two nodes in the u^ 
direction, where, we denote, for example, a state as (-,+), if it is odd under translation by 
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FIG. 2: M2 in eV, Lk = 0.5, L/2u; = 30, 27r, vr, 27r/3 



L, and even under reflection. 



B. Triple twist Mobius strip: Lk = 1.5 
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Shown in Fig. H] is an unsealed triple twist Mobius structure, constructed solely from the 
values of k and r of its centreline in [31]. The corresponding contour plot of M^ for the 
developed strip is shown at the top of Fig. [5l with the length scale arbitrary, along with 
three larger widths, including one at self-contact. The symmetry arguments developed above 
for Lk = 0.5 carry over for Lk = L5 with the replacement L -^ L/?). Thus K,,i],ri' have 
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FIG. 3: E',A', Lk = 0.5 

the same reflection symmetry as in Fig. [H with L -^ -^/3, and the same symmetry under 
the translation s ^ s + L/3. These induce corresponding parity eigenstates (I7l),(l9]) with 
L replaced with L/3. These are the so-called basically periodic solutions l4], which allow 
the domain to be reduced to L/6 for numerical integration using the FD method. From the 
point of view of Bloch's theorem, x(u^ + L/3, m^) = exp(i/cL/3)x(M^, — m^), with a Bravais 
lattice vector of magnitude L/3. To satsify the periodic boundary condition kL must be an 
integer multiple of vr. Here we only show solutions for kL/?) = , it. 



As with Mathieu's equation, solutions of higher period exist [l4|, so long as the periodic 
boundary condition is satisfied. For example, to show solutions symmetric or anti-symmetric 
under translation x(m^,m^) = ±x(m^ + L, —u^) are induced by the underlying symmetry of 
K, 7], 7]' , note that under s -^ s + L, u^ -^ u^ + L and that M^ is invariant (?7, rj' and 
u^ all change sign). Similarly, under the reflection, s-^L — s^u^^L — u^ and M^ is 
invariant since r] is even, but rj' and v? change sign under translation. From the point of 
view of Bloch's theorem, this corresponds to a Bravais lattice vector of magnitude L. For 
these higher periods fewer nodes are spread over the same length L, so more states of lower 
energies are observed. Note that period L/3 implies period L but not vice versa. Shown in 
Fig. [6] (left) is the dimensionless energy E' for Lk = 1.5. The large w energy behaviour does 
not now decrease faster than w"^. The (-,-) and (-,+) states also show a degeneracy at large 
w. 




FIG. 4: Lk = 1.5, unsealed L/2w = 9.36 structure with M"^ and generators shown. 



C. Curvature trapping of states 



Much attention has been given to the degree of locaUsation of states of particles con- 
strained to move in structures with curvature p, [l5|], with Costa [8] giving the example of 
one bound state for the bookcover surface with a zero transmission coefficient in the case 
of an infinitely sharp bend. One measure used for distinguishing between localised and 
extended states is the inverse participation ratio |7|, [l6|], A = J \^p\'^dV/{J \^p\'^dV)'^. The 
non-degenerate fiat ground state x = cos{ttu'^/2w) has dimensionless A' = A{2Lw) = 1.5, 
whereas the fiat higher excited states have A' = 2.25. 

Shown in Fig. [3] (right) is A' for Lk = 0.5, showing a general trend of increase in localisa- 
tion for increasing width. Comparison with the fiat state values is perhaps not so meaningful 
for higher excited states as they all have the same inverse participation ratio. The figure 
shows A', following the lowest energy eigenstate for each parity with increasing width. 

Shown in Fig. [6] (right) is the corresponding A' = A{2Lw/3) for Lk = 1.5. The fiat values 
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FIG. 5: M^ in eV, Lk = 1.5, L/2w = 9.36,4.72,4.49,4.22. 



are the same as for Lk = 0.5. There is a marked increase of the inverse participation ratio 
with width, except for the (-,+) state. There is a dip in the locahsation of the ground state 
before it also sharply increases with width. 

The corresponding lowest energy wave functions are shown in Fig. [7] for the highest 
width Mobius structures, showing confinement of the wave function to the high curvature 
regions corresponding to the lowest plot in each of Figs. [21 [51 The difference in topology 
is that Lk = L5 structures have no creases, prevented by self-contact, giving degenerate, 
disconnected, wave functions concentrated at the singularities in M^ at the largest width. 
The Lk = 0.5 structures, by contrast, show locahsation to the creases formed at higher 
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FIG. 6: E',A', Lk = 1.5 
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widths, allowing the wave function to be connected across the whole domain. 
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IV. CONCLUSIONS 

The eigenvalues and eigenstates are calculated for a quantum mechanical particle confined 
to geometric Mobius strips of arbitrary width, with Lk = 0.5, 1.5, the shape calculated by 
rigorously minimising the Wunderlich bending energy functional. The eigenstates are found 
to be localised to the potential wells obtained when increasing the width of the Mobius 
structures, leading to increased localisation to the regions of highest curvature. Our geo- 
metric formulation could be used to study transport properties of Mobius strip nanoribbons 



and molecules used in nanoscale devices 13| 
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